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ABSTRACT 

We discuss outer gap closure mechanism in the trans-field direction with the 
magnetic pair-creation process near the stellar surface. The gap closure by the 
magnetic pair-creation is possible if some fraction of the pairs are produced with 
an outgoing momentum. By assuming that multiple magnetic field will affect the 
local field near the stellar surface, we show a specific magnetic field geometry near 
the stellar surface resulting in the outfiow of the pairs. Together with the fact that 
the electric field is weak below null charge surface, the characteristic curvature 
photon energy emitted by incoming particles, which were accelerated in the outer 
gap, decreases drastically to ~ lOOMeV near the stellar surface. We estimate the 
height measured from the last-open field line, above which 100 MeV photons is 
converted into pairs by the magnetic pair-creation. We also show the resultant 
multiplicity due to the magnetic pair-creation process could acquire Me± ~ 10^ — 
10^. In this model the fractional outer gap size is proportional to P~^/^. The 
predicted gamma-ray luminosity (L^) and the characteristic curvature photon 
energy (Ec) emitted from the outer gap are proportional to B^P~^^^ and B^^^P~^ 
respectively. This model also predicts that and are related to the spin down 
power [Lsd) or the spin down age of pulsars (r) as oc L^^f or oc r"^/'', and 
Ec oc L^^'' or Ec oc r~^/^ respectively. 

Subject headings: pulsars: general- radiation mechanisms:non-thermal- -gamma 
rays:theory-magnetic field 
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1. Introduction 

The mechanism of particle acceleration and high-energy emission processes in the 
pulsar magnetospheres are one of the unresolved physics of the pulsar activities. The 
particle acceleration process and resultant high-energy 7-ray emission process have been 
discussed with the polar cap model (Ruderman & Sutherland 1975; Daugherty & Harding 
1982), the slot gap model (Arons 1981; Harding, Usov & Muslimov 2005; Harding et al. 
2008) and the outer gap model (Cheng, Ho & Ruderman 1986a,b; Hirotani 2008; Takata & 
Chang 2009). The polar cap model assumes the emission site close to the stellar surface 
above the polar cap, and the slot gap and the outer gap models assume the emission site 
in the outer magnetosphere. The different acceleration models have predicted the different 
properties of the 7-ray emissions from the pulsar magnetospheres. 

The observations of the pulsar emitting electromagnetic radiation in the high-energy 
7-ray bands have been facilitated by recent space and ground based telescopes. In 
particular, the Fermi 7-ray telescope has measured the 7-ray emissions from ~ 46 pulsars 
(Abdo et al. 2010, 2009a,b), including 21 radio-loud, 17 radio-quiet and 8 millisecond 
pulsars. In addition, AGILE (Astro-rivelatore Gamma a Immagini LEggero) has also 
reported the detection of 7-ray emissions from 4 new pulsars with 4 candidates (Pellizzoni 
et al. 2009). In more higher energy regime, MAGIC (Major Atmospheric Gamma Imaging 
Cherenkov) telescope has detected for the first time pulsed gamma-ray radiation above 
25 GeV from the Crab pulsar (Aliu et al. 2008). These observations will be useful to 
discriminate between the emission models. For example, the Fermi telescope has measured 
the spectral properties above 10 GeV with a better sensitivity than EGRET. It was found 
that the spectral shape of 7-ray emissions from the Vela pulsar is well fitted with a power 
low (photon index F ~1.5) plus exponential cut-off {Ecut ~ 3 GeV) model. The discovered 
exponential cut-off feature predicts that the emissions from the outer magnetosphere (Abdo 



-4- 



et al. 2009c) is more favored than the polar cap region (Daugherty & Harding 1996), which 
predicts a super exponential cut-off with the magnetic pair-creation. Furthermore, the 
detection of the radiation above 25 GeV bands associated with the Crab pulsar has also 
predicted the high-energy emission in the outer magnetosphere (Aliu et al. 2008). 

The pulse profiles observed by the Fermi telescope allow us to study the site of the 
7-ray emissions in the pulsar magnetosphere. Venter et al. (2009) fitted the pulse profiles 
of the 8 millisecond pulsars with the geometries predicted by the different emission models. 
They showed that most of the pulse profiles can be best fit with the outer gap (Takata et 
al. 2007; Takata & Chang 2007; Tang et al. 2008) or the two pole caustic (Dyks & Rudak 
2003; Dyks et al. 2004) geometries, which have a slab like geometry along the last-open 
field lines. However, they also found that the pulse profiles of two out of eight millisecond 
pulsars cannot be fitted by either the geometries with the outer gap or the caustic models. 
They proposed a pair-starved polar cap model, in which the multiplicity of the pairs is not 
high enough to completely screen the electric field above the polar cap, and the particles 
are continuously accelerated up to high altitude over full open field line region. 

The increase of the 7-ray pulsars allows us to perform a detail statistical study of the 
7-ray pulsars. In particular, the Fermi 7-ray pulsars including millisecond pulsars will 
reveal the relation between the 7-ray luminosity (L^) and the spin down power [Lsd), for 
which (X L^^ with f3 ~ 0.5 — 0.6 was predicted by EGRET measurements (Thompson 
2004). Also, the Fermi 7-ray pulsars will enable us to discuss the general trend of the 
relation among the spectral properties of the 7-ray emissions (e.g. the cut-off energy 
and photon index) and the pulsar parameters (e.g. rotation period and surface magnetic 
field). Together with the observed pulse profiles and the spectra, these general properties 
of the 7-ray emissions will discriminate between the 7-ray emission models in the pulsar 
magnetospheres. 
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In this paper, we discuss the 7-ray emissions from the outer gap accelerator. We 
propose a new outer gap closure mechanism by the magnetic pair-creation process near the 
stellar surface. The pairs produced by the magnetic pair creation will be able to close the 
gap if the sufficiently strong surface multiple field exists and affects the dipole field near the 
surface. In section O we first summarize results of gap closure process by photon-photon 
pair-creation process, and then we discuss our new gap closure mechanism by the magnetic 
pair-creation process. In section [3l we describe the model predictions of the properties of 
the 7-ray emissions. In section |U we will compare the model predictions with the results of 
the Fermi observations. We also discuss applicability of our model. In section O we will 
summarize our gap closure model and predictions for the outer gap accelerator. 

2. Gap closure mechanism 

2.1. Photon-photon pair-creation process 

The outer gap accelerator model was proposed by Cheng, Ho and Ruderman (1986a,b), 
who argued that a large global current fiow through the outer magnetosphere causes a charge 
depletion from the Goldreich- Julian charge density, which is defined by pcj = —Q ■ B/2'itc 
(Goldreich & Julian 1969) with Q being the vector of the rotation axis and B the magnetic 
field. In the charge depletion region, the non-corotational electric field along the magnetic 
field accelerates the charged particles, which result in the high energy 7-ray emissions. This 
non-corotational electric field could be screened out by the discharge of the copious electron 
and positron pairs produced by the pair-creation process of the 7-ray photons. The outer 
gap would be completely screened out in the trans-field direction, where the non-corotational 
electric field perpendicular to the magnetic field in the poloidal plane is equal to zero. 
This condition implies that the total potential (corotational -|- non-corotational) field is 
continuously connected to the corotational field outside the outer gap. 
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Zhang & Cheng (1997) discussed the gap closure mechanism by the photon-photon 
pair-creation process between the high-energy 7-rays emitted in the gap and the X-rays 
coming from the stellar surface. They estimated the typical gap thickness from the 
pair-creation condition E^Ex ~ 2(meC^)^, where E^ is the energy of the emitted 7-ray 
photons in the outer gap and Ex is the energy of the soft-photons from the stellar surface, 
and they obtained the fractional gap thickness as fp = h±{Ric/2)/ Ric ~ d.^P'^^^'^^B^^^'^ , 
where Ric is the light cylinder radius, h±{Ric/2) is the gap thickness in poloidal plane at 
r = Ric/2 and -8^,12 is the global stellar magnetic field in units of 10^^ Gauss. 

The outer gap closure mechanism with the photon-photon pair-creation process in the 
trans-field direction have also been discussed by solving the electrodynamics in the outer 
gap with 2-dimensional and 3-dimensional geometry (Takata, Shibata and Hirotani 2004; 
Hirotani 2006a,b; Hirotani 2008). For example, Hirotani (2006a) demonstrated that the 
outer gap for the young pulsar, the Crab pulsar, is almost screened out in the trans-field 
direction at the fractional gap thickness of / ~ 0.2 (figure 6 in Hirotani 2006a). For mature 
pulsars such like Gemiga, on the other hand, Takata and Chang (2009) argued that the 
photon-photon pair-creation process will be insufficient in the outer magnetosphere and the 
outer gap could occupy entire region between the last-open field lines and the critical field 
lines that have the null charge point at the light cylinder. 

2.2. New gap closure mechanism; magnetic pair-creation process near the 

stellar surface 

In this paper, we propose a possible gap closure mechanism, in which the magnetic 
par-creation process near the stellar surface supplies the electron and positron pairs to close 
the outer gap in the trans-field direction. It has been proposed that the magnetic field 
near the stellar surface is enhanced by the strong multiple magnetic field, although the 



- 7- 



global magnetic field is well described by the dipole field (Blandford et al. 1983; Romani 
1990; Ruderman 1991). The neutron star magnetic field will be produced by a current 
flowing the crust, which has a thickness of (5r ~ 1 — 3 x 10^ cm <^ Rs- It was suggested 
that the magnetic structure near the stellar surface is super position of clumps resolved 
into multiples of characteristic of the order of Rg/Sr (Arons 1993; Zhang & Cheng 2003). 
If the clumps cover whole stellar surface, the strength of the stellar magnetic field is order 
of Bg ~ {Rs/6r)"'Bd, where Bd is the strength of dipole magnetic field determined by the 
observed rotation period P and the period derivative P, and n = 1 and 2 represent coherent 
and incoherent superposition of the magnetic momentum of the clump, respectively. This 
model indicates the strength of the stellar magnetic field can take easily 5^ ~ 10 — lOO-B^. 
Therefore, even in the millisecond pulsar, the magnetic field near the star surface will be 
close to Bg ~ 10^^ Gauss, and the magnetic pair-creation process will take place close to 
the stellar surface. 

We will argue that the 7-rays emitted near the stellar surface is converted into the 
pairs via the magnetic pair-creation process above the height (/i_L,m) measured from the 
last-open field line. In equation (|6]), we will estimate the height, h±^rn, above which the 
pair-creation process takes place. The important conditions for closing the outer gap by 
the pairs produced by the magnetic pair-creation process near the stellar surface are as 
follows; (1) the gap was not closed by the photon-photon pair-creation process below the 
height h± m (c.f. section . (2) the local magnetic field lines near the stellar surface is 
bending away from the last-open field lines (see below and Figure [2]) due to the strong 
multiple fields, and (3) some of magnetic pairs migrate into outer magnetosphere. The 
second condition is required to produce outflows of the magnetic pairs. 

Figure [1] and [2] represent the schematic view of the outer magnetosphere and of the 
magnetic structure near the stellar surface, which is favored in this study, respectively. In 
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the outer magnetosphere, the photon-photon pair-production process create pairs, which 
are separated by the accelerating electric field. Inward propagating particles will emit 
the 7-ray photons toward the strong magnetic field region near the stellar surface. As 
demonstrated by the electrodynamic study (e.g. Hirotani 2006a), the electric field below 
null charge surface are significantly reduced by the pairs and arises with a very weak field. 
Below the null charge surface, therefore, the curvature energy loss will not be able to be 
compensated by the acceleration of the electric field in the gap (c.f. section IT3|) . and the 
incoming particles loose their energy by the curvature radiation. It is interesting to note 
that there is a minimum energy of the curvature photons, which does not depend on any 
pulsar parameters and the curvature radius of the local magnetic field. Assuming that the 
curvature loss dominates the energy gain due to the acceleration by the electric field, the 
evolution of the Lorentz factor may be described as 



where s is the curvature radius. Close to the null charge surface, the incoming particles 
lose their energy with a time scale smaller than the time scale of travelling to the stellar 
surface, because the Lorentz factor is high enough. As decreasing the Lorentz factor of the 
particles, the curvature energy loss time scale becomes comparable to the crossing time 
scale of (it ~ s/c. In such a case, equation of motion ([1]) implies 1/7^ — 1/7q ~ 4e^/3meC^ 
where 70 is the initial Lorentz factor. If 7 << 70, the typical Lorentz factor of the particles 
below the null charge surface becomes 



where S7 is the curvature radius in units of 10^ cm. 

The Lorentz factor ([2]) gives the minimum energy of the curvature photons, which does 
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(2) 
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not depend on any pulsar parameters, that 

3 h-f^c 9mec2 

£^mm, ~ 7 ~ — 77 MeV, (3) 

4 s 8a/ 

where «/ is the fine structure constant. We expect that the incoming electrons will mainly 
emit the curvature photons with an energy Emm ~ 100 MeV between the stellar surface and 
the null charge surface. Some of 100 MeV photons emitted below the null charge surface 
pass through vicinity of the stellar surface. Applying the static dipole field geometry and 
P = 0.1, for example, the trajectory of the 100 MeV photons emitted below r ~ 10^ cm on 
the last-open filed line will cross the magnetic pole below r ~ 2 x 10^ cm. Therefore, it is 
expected almost all 100 MeV photons emitted below r ~ 10^ cm will be converted into pairs 
with a strong magnetic field of 5 ~ 10^^ — 10^'^ Gauss. We find that an incoming particle 
emits about ~ 2 x lO^s^'^^'^ photons before reaching the stellar surface below r ~ 10^ cm . 

For the millisecond pulsars, the 100 MeV curvature photons could be converted into the 
pairs with a strong multiple magnetic field, which could acquire ~ 10^^ Gauss, at very 
close to the stellar surface. With the static dipole field geometry, the curvature photons 
emitted blow r ~ 3 x 10^ cm is propagating toward the magnetic pole for the typical 
rotation period of P = 4 ms. An incoming particle will emit A^ ~ lO^Sg of ~ 100 MeV 
photons below r ~ 3 x 10^ cm 

For the canonical pulsar, the pairs created by the magnetic pair-creation process 
of 100 MeV photons can also emit the soft 7-rays via the synchrotron radiation, 
which may further generate new pairs. With the typical Lorentz factor of F ~ 100, 
one pair emits about 1~10 synchrotron photons with a typical energy of Egyn ~ 
5(F/100)2(sin^j0.05)(5/5 ■ 10" G) MeV. This ~5 MeV synchrotron photons also could 
be converted into pairs with a strong local magnetic field 10^^ ~ 10^'^ Gauss near the stellar 
surface, where the magnetic field will be enhanced by the multiple magnetic field. This 
implies the resultant multiplicity of an incoming particles, which was accelerated in the 
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outer gap, could be Me± ~ lO^s^^''^ for the canonical pulsars. 

The mean free path of the magnetic pair-creation process of the photon with the energy 
is described as (Erber 1966;Ruderman & Surtherland 1975) 

4.4 h Bg f 4:\ 



lm = - ^exp — , (4) 

afUiecB^ \oX/ 

where x = E^B^/2m.f,c^Bq, Bg = 4.4 x 10^^ Gauss and B± = BsinOa with 9a is the 
angle between the direction of the propagating for the photon and the magnetic field. 
The exponential dependency of the mean-free path on the photon energy E^y and the 
angle 6a implies that most of emitted photons will be converted into pairs if the condition 
that E^B ^/ {2mf,c^ Bg) ~ x is satisfied. Using the magnetic pair-creation condition, we 
can estimate height measured from the last open field line, above which the magnetic 
pair-creation process of the photons with Emin ~ 100 MeV becomes important process. 
The collision angle 6a is approximately described as 



sin6la ~ ^ ~ (5) 



where i ~ \/2h±s is the propagating distance of a photon from the emission point, and 
/i^ is height measured from the last-open field line. The magnetic pair-creation condition 
implies that the curvature photons with an energy ~ 100 MeV will be converted into the 
pairs above the height 

h±,m{.Ri) ~ 10^X-i^m%-57 cm, (6) 

where Ri expresses the critical radial distance below which the magnetic pair-creation 
process becomes to be important and it will be -Rj ~ 2 — 3i?s for the canonical pulsars 
and Ri ~ Rs for the millisecond pulsars. In addition, x~i = x/0-1, -Bm,i2 is the strength 
of the magnetic field at the pair-creation position in units of 10^^ Gauss, and we used 
Emin = 100 MeV. Rescaling the thickness at the stellar radius using the magnetic flux 
conservation, we obtain h^^rn{Rs) ~ {Rs/ RiY^'^h±^ra{Ri) with Rs being the stellar radius. 
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As we have discussed above, the incoming particles will create pairs above the height 
h±,m, which is described by equation dH]), via the magnetic pair-creation process near the 
stellar surface, and the multiplicity could acquire M^± ~ 10^ — 10^ for the canonical pulsars 
and ~ 10^ for the millisecond pulsars. Although it is expected most of the created pairs 
have inward momentum and migrate toward the star, it may be possible that small fraction 
of the created pairs are produced with an outgoing momentum and migrate into outer 
magnetosphere, as discussed below. In fact, about ~ 10 pairs out of Me± ~ 10^ — 10^ will 
be enough to close the outer gap accelerator in the trans-field direction everywhere in the 
outer magnetosphere. 

The magnetic pairs would be produced with an outgoing momentum if the magnetic 
structure near the stellar surface is affected by the strong multiple field, as illustrated in 
Figure [21 which shows how some pairs can acquire the outgoing momentum due to the 
geometry of the local magnetic field lines. In fact, it is required that the magnetic field 
(thick dashed-lines) on the pair-creation points are bending away from the last-open field 
line due to a strong multiple field (solid-dashed line). The solid curved lines represent 
the global field which are not affected by the local magnetic field, and thin dashed-lines 
represent the geometry without the multiple field. The 100 MeV curvature photons (solid 
arrows) emitted by the incoming particles, which were accelerated inside the outer gap, are 
converted into the pairs via the magnetic pair-creation above height, h±^rn described by 
equation ([6]). For the canonical pulsars, the synchrotron radiation of the created pairs will 
produces ~5 MeV photons (dashed arrows), which could be converted into new pairs by a 
strong magnetic field near the stellar surface. If the direction of the magnetic field at the 
pair-creation positions is bending away from the last-open field line, it is possible that the 
collision angle (denoted as a in the Figure [2]) between the photons and the magnetic field 
line becomes larger than a > 90°, implying the created pair has the outgoing momentum. 
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We note that (1) because the position of the pair-creation point is determined by the 
value i?sina, and (2) because the coUision angle changes from a smaller value to a larger 
value as the photon propagate toward the stellar surface, the position with the collision 
angle smaller than 90° is preferentially chosen as the pair-creation point, although sin a 
gives the same values with, for example, a = 80° and 110°. However, the magnetic field also 
increases as the photons propagate toward the stellar surface, implying the strength of the 
magnetic field is bigger for the position with a = 110° than that with a = 80°. As a result, 
it will be possible that some photons do not have enough energy to produce the pairs with 
the magnetic field at the position with a = 80°, but do to produce the pairs at the position 
with a = 110°. On these ground, some created pairs will be able to have the outgoing 
momentum, if the local magnetic field is bending away from the last-open field line. 

With the magnetic field structure bending away from the last-open field line, the 
incoming particles from the outer gap could emit the 100 MeV curvature photons toward 
lower altitude, and therefore could make pairs at lower altitude, implying the magnetic 
pair-creation process could occur inside the gap if the inner boundary is located close to 
the stellar surface. The created pairs will be discharged by the small electric field along 
the magnetic field and this discharge will partially screen the electric field. Because the 
created pairs lose their perpendicular momentum within very short distance, the Lorentz 
factor after loosing the perpendicular momentum is F ~ 1/ sin 6*0 ~ 10(sin ^q/O.I)"^. The 
magnetic pairs will screen the electric field near the inner boundary so that the potential 
drop between the inner boundary and the point, below which the magnetic pair-creation 
process is occurred, is reduced to ~ 5 x 10^ Volts. 

If all field lines in the polar cap region are bending toward the last-open field line, the 
outgoing magnetic pairs are not expected because the collision angle is always a < 90°. 
In such a case, all pairs created by the magnetic pair-creation process will have an 



- 13 - 



inward momentum, implying the photon-photon pair-creation process will be only possible 
mechanism to close the outer gap. 

3. Outer gap closed by magnetic pair-creation process 

Our gap closure process is summarized as follows; 

1. the incoming particles emit ~100 MeV photons below the null charge surface, 

2. the ~100MeV photons emitted toward the stellar surface will produce the pairs via 
the magnetic pair-creation process (Me± ~ 10^ — 10^) above the height h±^rn estimated 
by equation and 

3. if some fraction of the created pairs will be produced with an outgoing momentum 
due to the geometry of the local magnetic field, then the outgoing pairs close the 
outer gap in the trans-field direction everywhere in the outer magnetosphere. 

We will discuss the applicability of the present model in section 14.21 

Applying equation as the typical thickness of the gap rescaled at the radial distance 
r = Ri with the dipole field geometry, we define the fractional gap thickness on the stellar 
surface as 

fm ^ ^^^^ ~ 0.25K{x, B.^, s)P'_{\ (7) 

with 

K = x\B-^%s, I^^J (8) 

Applying the fractional gap thickness to the acceleration and curvature emission 
process beyond the null charge surface, we discuss the expected properties of the 7-ray 
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radiation from the outer gap closed by the magnetic pair-creation process. The typical 
magnitude of the electric field in the gap beyond the null charge surface is given by 

E„(i^, P) ~ 4^ ~ 8.3 X IQ^K'^B^^.^PZI Volt /cm. (9) 

which can accelerate the electrons and positrons up to 

1{K, B,) ~ ~ 2.1 X 10^/^^251/4 ^ (10) 

where we used s = Ric in the outer magnetosphere. We find that the maximum Lorentz 
factor does not depend on the rotational period. The typical energy of the curvature 
radiation by the accelerated particles in the outer magnetosphere is 

E,{K, B,,P) = ^^^ 0.55K'/'Bl{\Pl^ GeV. (11) 

Because the efficiency of the emission is significantly reduced above the energy Ec, we 
expect that the 7-ray spectrum has the cut-off energy of Ec- 

The total current flowing in the gap is order of Igap ~ fmlcj, where Iqj = irBdRs/PRic 
is the Goldreich-Julain current. Using the total potential drop in the gap of Vgap ~ fm^a, 
we can estimate the 7-ray luminosity as 

L^{K, 5,, P) ~ IgapVgap ~ 2 X IO^^K^BI^^P:^/^ erg/s. (12) 

4. Application and discussion 

4.1. Comparison with Fermi observations 

The present model predicts the cut-off energy and the luminosity of the 7-ray emissions 
are related with the stellar magnetic fleld and the rotation period as Ec oc B^J^P-^ of 
equation ( ITT]) and oc BjP~^^'^ of equation (IT^ . respectively. We plots the our predictions 
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with the observations on the cut-ofF energy in Figure [3] and on the luminosity in Figure HI 
In the figures, the circle and triangle symbols represent the radio-selected and 7-ray selected 
7-ray pulsars, and the squares are the millisecond pulsars. We can see that the slope of 
the model prediction is consistent with the observations for both canonical pulsars and 
millisecond pulsars. 

In the present outer gap model, the properties of the 7-ray emissions depend on the 
local parameters K{x,Bm,s) = X-iB^\2^7iRs/ Ri)^^"^ , which is determined by the local 
magnetic structure. In fact, we adopted K = 2 for the canonical 7-ray pulsars and K = 15 
for the millisecond pulsars in Figures |3] and IH The strength of the multiple field at a 
distance SR from the stellar surface may be expressed as Bm ~ Bs[{5r + SR) /Sr]~^'^~^^\ 
where Bg is the strength of the multiple field at the stellar surface and 5r ~ 1 — 3 x 10^ cm 
is the thickness of the crust, the index m is the multiplicity. If we consider the localized 
dipole field (m = 2), B^ ~ Bs[{Sr + 6R)/6r]~^, with B^ ~ 10^^ Gauss, this localized dipole 
field becomes same order of magnitude with the global dipole field at about one stellar 
radius from the stellar surface, because Bm{r — Rg = Rg) ~ Bs{6r / RsY ~ 10~^-Bs with 
5r ~ 2 — 3 X 10^ cm and because the global dipole field becomes B^jl? at r — Rg = Rg. 
If the radial distance below which the magnetic pair-creation process takes place is 
Ri = 2Rg, the strength of the magnetic field at Rt is Bm. ~ 10^^ Gauss and the curvature 
radius of the multiple field will be s ~ 10^ cm. This implies the local parameter of 

^^2-57 ~ 3.5, which explains K ~ 2 in Figures |3] and HI For the millisecond 
pulsars, it is expected that the magnetic par-creation of the photons with the energy 
Emin ~ 100 MeV is possible only at very close to the stellar surface (r ~ Rg), because the 
field strength is about three or four order of magnitude smaller than that of the canonical 
pulsars. If we adopt i?s ~ 3 x 10^'^ Gauss as the strength of the multiple field on the stellar 
surface and s ~ 10^ cm as the curvature radius of the multiple field, we obtain K ~ 10. 
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As Figures |3] and H] show, the present model predicts that the local parameter 
K{xi Bmi s) discriminates the canonical and the millisecond pulsars as the different 
populations on the plots of the cut-off energy E^. versus B^J^P~^ and of the luminosity 
versus i?JP~^/^. From equations ( 1111) and (fT2l) . on the other hand, one can see that Lry/El 
does not include the local factor K{x, Bm, s). More strictly speaking, the quantities L^/E^ 
carries away the dependency on the gap thickness /, which depends on the local parameters 
K, because oc f^Lsd and Ec oc f^^'^Va^^s~^^'^. We then obtain the relation that 

^ ~ 6.6 X 10''bI%P:I erg/s GeVl (13) 

Figure [5] compares the model prediction and the observations for each pulsar. Figure |5] 
shows that both canonical pulsars and millisecond pulsars are consistent with the line of 
L^l E^ oc B^^^P-i- This predicts that although there is a gap between the two populations, 
where no 7-ray pulsars are plotted in Figure O in fact two populations will be continuously 

1 /2 1 

connected. It is expected that more 7-ray pulsars having smaller BJ^2P-i than the present 
canonical 7-ray pulsars will be discovered by Fermi 7-ray telescope. Those pulsars will be 
expected to distribute around solid line in Figure El and two populations, i.e. the canonical 
and millisecond pulsars, will be continuously connected in the plot. 

Recasting the cut-off energy E^. given by equation (fTTl) and the 7-ray luminosity 
given by equation (fT2l) as a function of the spin down age r = P/2P, we obtain 

~ 7K'''B-li\^'" GeV (14) 

and 

L, ~ IQ''K'B~,%\;'/' erg/s, (15) 

respectively, where is the spin down age in units of 1 kyrs and ^3^,34 is the spin down 
power in units of 10'^^ erg/s. Figures [6] and [7] plot the cut-off energy and the 7-ray luminosity 
as a function of the spin down age, respectively. The lines in Figure [6] and Figure [7] are 
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results for K = 2 and the typical magnetic field of < -8^,12 >= 3 for the canonical pulsars 
(solid lines), and for K = 15 and < -8^,12 >= 3 x 10""^ for the millisecond pulsars (dashed 
lines). We can see that the model predictions are consistent with the observations for both 
populations. The present model predicts that the two populations will be separated in the 
plots of the emission properties {Ec, L-y) as a function of the spin down age. 

If we recast the cut-off energy Ec given by equation f lTT]) and the 7-ray luminosity 
given by equation ( 1T2|) as a function of the spin down power Lgd = (27r)^i?^i?^/6c^P^, we 
obtain 

Ec ~ 0.31K'/'B]%L% GeV (16) 

and 

~ 5 X IQ-^' K'b'J,\L%%, erg/s (17) 

respectively. One can see that the proportionality coefficients K^^'^B^^^ equation (fT6l) 
and K^B^^^2 equation f lT7|) have less dependency on the pulsar populations. For example, 
if we apply K = 2 and the typical field < -8^,12 >= 3 for the canonical pulsars and K = 15 
and < Bd i2 >= 3 x 10"^ for the millisecond pulsars, we obtain < K^/^B^^^ >~ 3.7 for 
the canonical pulsars and ~ 7.6 for the millisecond pulsars, and < K^B^J^2 >~ 18 the 
canonical pulsars and ~ 7.7 for the millisecond pulsars. We find that the proportionality 
coefficients between the two populations are different only about factor of two, implying it 
is difficult to discriminate between the two populations in the plot of the cut-off energy (or 
7-ray luminosity) versus the spin down power. Figure |8] and [9] plot the cut-off energy and 
the 7-ray luminosity as a function of the spin down power, respectively. The solid lines are 
model predictions with < K^/'^B]^^2 >~ 5 in Figure m and < K^B'^J^^ 

>~ 13 in Figure |9l 

Figures |8] and [9] show the predicted relations that Ec oc L^J^ and oc E'J^ are consistent 
with the observations for both canonical pulsars and millisecond pulsars. 
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4.2. Applicability to pulsars 

In this section, we describe applicability of our model to the pulsars. The present gap 
closure model invokes the condition that the thickness h^^^, above which the magnetic 
pair-creation becomes to be important, is less than the thickness, at which the gap is 
closed by the photon-photon pair-creation process. To discuss the applicability to the 
pulsars, for example, one may compare the present model with a self-consistent model 
investigated by Zhang & Cheng (1997), who discussed the gap thickness closed by the 
photon-photon pair-creation process. They estimate the fractional gap thickness, which is 
defined by the ratio between the typical gap thickness at r = Ric/2 and the light radius, 
as fp = h±{Ric/2) / Ric ~ 

our fractional gap thickness defined in 
equation ([7]) into one in the sense of the definition by Zhang & Chang (1997). Because the 
trans-field thickness of the magnetic fiux tube is approximately proportional to h± oc r^/^, 
the fractional gap thickness is described by fm,i = h^^rn{Ric/'^) / Ric ~ 2"^/^/^, where fm is 
defined by equation ([7]). 

We compare the predicted thicknesses of fp and fm^i for 36 canonical 7-ray pulsars in 
Table 1 and for 8 millisecond 7-ray pulsars in Table 2. We find that the pulsars with a 
larger spin down luminosity, such like PSR J0534-I-2200 (the Crab pulsar), has fp« fm,i- 
As Table 1 and Table 2 show, on the other hand, the mature pulsars, such like PSR 
J0633-I-1746 (the Geminga pulsar), and the millisecond pulsars, who show a smaller spin 
down power, indicate fp» /m,,i- 

The present gap closure process with the magnetic pair-creation process will be 
applicable for the pulsars whose have fm,i ^ fp, that is, the predicted gap thickness fm,i 
is comparable with or slightly less than fp of the photon-photon pair-creation process. 
The condition that fp « fm,i implies that the efficient photon-photon pair-creation will 
definitely close the gap before the gap reaches the thickness fm,i ■ For the case fp» fm,i, 
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on the other hand, the outer gap must be thicker than /„, but the gap may be closed by 
magnetic pair-creation process with the thickness at which the photon-photon pair-creation 
is initiated in the gap. Equating the fractional gap thickness fp and /m,i, it may be 
suggested that the switching of the gap closure mechanism takes palace at 

L,, ~ 2.5 X 1036(/^/2)-^6«/3ii?,;?f ^ erg/s, (18) 

indicating the photon-photon pair-creation closes the gap for the pulsars with Lgd ^ 
]^g36 erg/s, while the magnetic pair-creation may close the gap for Lsd ^ 10^^ erg/s. 



4.3. Validity of assumption 

In the present model, we assume that the energy loss due to the curvature radiation 
below the null charge surface is not compensated by the acceleration of the electric field in 
the gap. For simplicity, we assume that the electric field decreases quadratically along the 
field lines below the null charge surface toward the inner boundary; 

where Rin and Rn are the radial distance to the inner boundary and the null charge point, 
respectively, and E\\fi is the electric field at the null charge point, which is characterized by 
equation ([9]). If the inner boundary is located ~ 1 stellar radius from the stellar surface, 
the acceleration field is reduced to E\\{r ~ R^) ~ {Rin/ Rn)'^E\\ Q ~ 10~^£'||^o near the inner 
boundary. Applying the typical electric field described by equation at the null charge 
point, a particle gains energy by the electric field with a rate of 

eE\\{r ~ Rin)c ~ iOK^Ba^uPli erg/s. (20) 

On the other hand, the particles having the Lorentz factor expressed by equation (j2]) loose 
their energy via the curvature radiation with a rate of 

Pc = ^7'^ ~ 5 X lOh-^^' erg/s, (21) 
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indicating the energy loss by the curvature radiation dominates the energy gain by the 
acceleration of the electric field. In fact, (1) the electrodynamic model expects the electric 
field will decreases more rapidly below the null charge surface than the quadratic expression 
(Hirotani 2006a) or (2) the inner boundary will located far from the stellar surface for 
the small current case. Therefore we can safely assume that the energy loss due to the 
curvature radiation below the null charge surface is not compensated by the acceleration of 
the electric field. 



4.4. Comparison with the previous works 

We briefly summarize the difference in the predicted 7-ray luminosity between the 
present model and the model investigated by Zhang & Cheng (1997). Zhang& Cheng (1997) 
predicted the fractional gap thickness fp ~ 5.6P^^/^^i?^ of the outer gap closed by the 
photon-photon pair-creation process in the outer magnetosphere. With this model, the 
7-ray luminosity depends on the spin down power as oc B)!'^ . Therefore, the 7-ray 
luminosity is less dependent on the spin down power. Later, they developed their gap model 
by taking into account the effects of the inclination angle between the rotation axis and 
the magnetic axis (Zhang et al. 2004). They argued that the fractional gap thickness and 
the resultant 7-ray luminosity depend on the inclination angle. They assumed that when 
the fractional gap thickness goes to unity, then the 7-ray luminosity approaches to the spin 
down power. Therefore, the 7-ray luminosity bounded between L^^^^ and L<jrf depending 
on the inclination angle. They carried out a Monte Carlo simulating to calculate the 7-ray 
luminosity by assuming that the inclination angle is randomly distributed. As a result, 
they obtained the relation that oc with an index of /3 = 0.38 ~ 0.46. In the present 
model, on the other hand, the fractional gap thickness oc P^/^ has a less dependency 
on the rotational period compared with the model of Zhang & Cheng (1997), in which 
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fp oc p26/2i ];gguit, without introducing the effects of the inchnation angle, the present 

model predicts the relation of oc L^^ with more steep index of /3 = 5/8. 

5. Conclusion 

In this paper, we have studied the outer gap accelerator model closed by magnetic 
pair-creation process. We argued that below null charge surface, the curvature loss is 
not compensated by the acceleration due to the electric field in the gap. In such a 
case, the incoming particles, which were produced in the outer magnetosphere, will emit 
curvature photons with about Emin ~ rriec'^/af ~ 100 MeV. The 100 MeV curvature photos 
propagating toward the stellar surface will be converted into pairs by the pair-creation 
process with the strong local magnetic field near the stellar surface, where the multiple 
magnetic field affects to the global field lines. For the canonical pulsar, the synchrotron 
radiation of the created pairs produce ~5 MeV photon, which will be furthermore converted 
into pairs. As a result, multiplicity of an incoming particle could acquire Me± ~ 10^ — 10^. 
With the local field lines bending away from the last-open field line (such as illustrated in 
Figure El), the created pairs via the magnetic pair-creation process can have the outgoing 
momentum and migrate into the outer magnetosphere. If ~ 10 pairs out of Me± ~ 10^ — 10^ 
migrate outward to the outer magnetosphere, those pairs could close the outer gap 
accelerator. According to this scenario, the main results of this paper are as follows. The 

1/2 

fractional thickness of the outer gap becomes fm ~ 0.25KPJi , which has a less dependency 
on the rotational period compared with the outer gap model proposed by Zhang & Cheng 
(1997). With the present model, the spectral properties of the 7-ray emissions depend on 
the local parameter K ~ x'^iB^i2S7{Rs/ Ri)^^"^ , which is determined by the local magnetic 
structure near the star. We expect that the local parameter K takes a vale of K ~ 2 for the 
canonical pulsars and K 15 for the millisecond pulsars (Figures |3] and H]). The present 
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model predicts that the canonical pulsars and the millisecond pulsars are connected in 
the plots of L.y/ Ef. versus B^^P"^/^, in which the effect of the fractional gap thickness is 
carried away (Figure E]). The present model predicts that the cut-off energy [Ec] and the 
7-ray luminosity (L^) depend on the spin down age or the spin down power as Ec oc r~^/^ 
and oc r~^/^ or Ec oc L^J^ and oc L'J^ (Figures [MSI)- 

In addition to the cut-off energy and 7-ray luminosity, which have been discussed in 
this paper, the Fermi 7-ray telescope provides the photon index of 7-ray spectrum and 
the pulse profiles for each pulsar (Abdo et al. 2010). It must be important to discuss the 
photon index and the pulse profile with the acceleration model, because they will contain 
information of the electric structure (e.g. the distribution of the electric field) in the 
acceleration region and the three-dimensional geometry of the emission region (e.g. Romani 
& Yadigaroglu 1995; Cheng, Ruderman & Zhang 2000; Spitkovsky 2006). However, a more 
detail model, which has to consider the electrodynamics in the gap and three-dimensional 
structure, is required to study the shape of the 7-ray spectra and the pulse profiles. 
Studying the emission properties with the electrodynamics in the present outer gap closure 
model will be done in the subsequent papers. 
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Name 


Pis) 


5d(10i2 G) 


L,,(1034 erg/s) 


fp 


fm,l] K = 2 


J0534+2200 


0.033 


3.8 


46100 


0.038 


0.10 


J1833-1034 


0.062 


3.5 


3370 


0.085 


0.14 


J0205+6449 


0.066 


3.6 


2700 


0.091 


0.14 


J2229+6114 


0.052 


2.0 


2250 


0.093 


0.13 


J1124-5916 


0.135 


10.0 


1190 


0.12 


0.21 


J1420-6048 


0.068 


2.4 


1000 


0.12 


0.15 


J0835-4510 


0.089 


3.4 


688 


0.14 


0.17 


J1813-1246 


0.048 


0.9 


626 


0.13 


0.12 


J1418-6058 


0.11 


4.4 


495 


0.16 


0.19 


J1952+3252 


0.040 


0.5 


374 


0.15 


0.11 


J1826-1256 


0.11 


3.7 


358 


0.17 


0.19 


J 1709-4429 


0.102 


3.0 


341 


0.17 


0.18 


J2021+3651 


0.104 


3.2 


338 


0.17 


0.18 


J1907+06 


0.107 


3.1 


284 


0.18 


0.18 


J1747-2958 


0.990 


2.5 


251 


0.19 


0.18 


J1048-5832 


0.124 


3.5 


201 


0.20 


0.20 


J1718-3825 


0.075 


1.0 


125 


0.22 


0.15 


J1459-60 


0.103 


1.6 


91.9 


0.25 


0.18 


J2238+59 


0.163 


4.1 


90.3 


0.26 


0.23 


J1028-5819 


0.091 


1.2 


83.2 


0.25 


0.17 


J1509-5850 


0.089 


0.9 


51.5 


0.29 


0.17 


J0007+7303 


0.316 


10.6 


45.2 


0.34 


0.31 


J1809-2332 


0.147 


2.2 


43 


0.32 


0.21 


J1958+2846 


0.29 


8.0 


35.8 


0.36 


0.30 
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Name 


P{s) 


BdilO'' G) 


LsdilO^^ erg/s) 


fp 


fm.,1] K = 2 


J2032+4127 


0.143 


1.7 


26.3 


0.37 


0.21 


J0248+6021 


0.217 


3.4 


21 


0.41 


0.26 


J0631+1036 


0.288 


5.5 


17.3 


0.45 


0.30 


J0742-2822 


0.167 


1.7 


14.3 


0.45 


0.23 


J1732-31 


0.197 


2.2 


13.6 


0.46 


0.25 


J0633+0632 


0.297 


4.8 


11.9 


0.50 


0.30 


J2021+4026 


0.265 


3.8 


11.6 


0.50 


0.29 


J2043+2740 


0.096 


0.3 


5.6 


0.55 


0.17 


J0659+1414 


0.385 


4.3 


3.8 


0.73 


0.35 


J0633+1746 


0.237 


1.6 


3.3 


0.71 


0.27 


J1057-5226 


0.197 


1.1 


3 


0.70 


0.25 


J1836+5925 


0.173 


0.5 


1.2 


0.92 


0.23 


J1741-2054 


0.414 


2.3 


0.9 


1.14 


0.36 


J0357+32 


0.444 


1.9 


0.5 


1.40 


0.37 



Table 1: Pulsar parameters, which were taken from Abdo et al. (2010), and the fractional 
gap thickness predicted by Zhang & Cheng (1997), /p, and by the present model, for 
38 canonical 7-ray pulsars. Here, the fractional gap thickness is defined by the ratio between 
the typical gap thickness at r = Rid 2 and the light radius. 
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Name 


Pirns) 


5^(10^ G) 


L,d(1034 erg/s) 


fp 


fmX, K=15 


J0218+4232 


2.3 


4.1 


24 


0.26 


0.20 


J0613-0200 


3.1 


1.8 


1.3 


0.60 


0.23 


J0751+1807 


3.5 


1.5 


0.6 


0.76 


0.25 


J1614-2230 


3.2 


1.2 


0.5 


0.78 


0.24 


J1744-1134 


4.1 


1.8 


0.4 


0.84 


0.27 


J2124-3358 


4.9 


2.4 


0.4 


0.89 


0.29 


J0030+0451 


4.9 


2.2 


0.3 


0.92 


0.29 


J0437-4715 


5.8 


2.9 


0.3 


0.98 


0.32 



Table 2: Pulsar parameters, which were taken from Abdo et al. (2010), and the fractional 
gap thickness predicted by Zhang & Cheng (1997), fp, and by the present model, /m,i, for 8 
millisecond 7-ray pulsars. Here, the fractional gap thickness is defined by the ratio between 
the typical gap thickness at r = Ric/2 and the light radius. 




Fig. 1. — Schematic view of the outer gap accelerator. In the outer magnetosphere, the 
photon-photon pair-creation process produce the pairs in the gap uniformly. The pairs 
created by the magnetic pair-creation process close the outer gap in the outer magnetosphere. 
The favorable magnetic structure near the stellar surface is depicted in Figure [2l 
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Global field 




Stellar surface 



Fig. 2. — Schematic view of the favorable magnetic structure near the stellar surface. 
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Millisecond . ■ Gamma-ray selected 
Model(K=15) Model {K=2) 




Fig. 3. — Plot of the cut-off energy as a function of B'^J^2P~i - The circle and triangle symbols 
are observations for the radio-selected and the 7-ray selected 7-ray pulsars, respectively. In 
addition, the square symbols represent the millisecond pulsar. The solid and dashed lines 
show the prediction of the model, E^. = 0.55K^^'^B^^^2P~i GeV of equation (ITT]) . The model 
results are for = 2 for the canonical pulsar (solid line) and K = 15 for the millisecond 
pulsars (dashed line). The observations are taken from Abdo et al. (2010). 
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Fig. 4. — Plot of the 7-ray luminosity as a function of i?J 12-^-1 • '^^^ symbols cor- 
respond to same case as Figure [31 The lines represent the model prediction, ~ 
2 X 10^^ K^B^ ^2^-1^^ erg/s of equation ( fT2l) . with K = 2 for the canonical pulsar (solid 
line) and i^' = 15 for the millisecond pulsars (dashed line). 




Fig. 5. — Plot of L^/E^ as a function of bJ^^P-i ■ The symbols correspond to 
same case as Figure [31 The line corresponds to the model prediction, L^/El ~ 6.6 x 
IQi^^B^/l^P-l erg/s GeV^ of equation 
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Fig. 6. — Plot of the cut-off energy as a fucntion of the spin down age. The sym- 
bols correspond to same case as Figure [31 The lines correspond to the model prediction, 
Ec ~ TJ'^^/^i?^ ^[2^X3 GeV of equation ([14]), with K = 2 and the typical magnetic field 
< -Brf,i2 >= 3 for the canonical pulsars (solid line), and K = 15 and < Bd^i2 >= 3 x 10~^ for 
the millisecond pulsars (dashed line). 
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Age (kyrs) 

Fig. 7. — Plot of the 7-ray luminosity as a function of the spin down age. The symbols 
correspond to same case as Figure [31 The lines correspond to the model prediction, ~ 
1^"^^ B ^\!^ T-j^ erg/s of equation ( fTSl) . with K = 2 and < -8^,12 >= 3 for the canonical 
pulsars (solid line), and K = lb and < 5^,12 >= 3 x 10""^ for the millisecond pulsars (dashed 
line) . 




Model (<k3'2b|S =5>) 

10"' V" ' \ ' ' ' ' ' ' 3 ' ' 4 ^5 

10 10 10 10 10 10 10 

Lgtj ; Spin down luminosity {10^^ erg/s) 

Fig. 8. — Plot of the cut-off energy as a fucntion of the spin down power. The sym- 
bols correspond to same case as Figure [31 The line is the model prediction, ~ 
0.22/^3/251/4^1/4^ GeV of equation with < K'^''^B]J^^^ >= 5. 
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Fig. 9. — Plot of the 7-ray luminosity as a function of the spin down power. The 
symbols correspond to same case as Figure [3l The line is the model prediction, ~ 
2 X 1032i\:3s^/4 ^5/8^ g^g/g equation (JIZD, with < K^B^^^^^ >= 13. 



